'Tilting' the Universe with the Landscape Multiverse: The 'Dark' Flow 
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The theory for the selection of the initial state of the universe from the landscape multiverse predicts super- 
horizon inhomogeneities induced by nonlocal entanglement of our Hubble volume with modes and domains 
beyond the horizon. Here we show these naturally give rise to a bulk flow with correlation length of order hori- 
zon size. The modification to the gravitational potential has a characteristic scale Lj ~ lO-'//^ ' , and it originates 
from the preinflationary remnants of the landscape. The 'tilt' in the potential induces power to the lowest CMB 
multipoles, with the dominant contribution being the dipole and next, the quadrupole. The induced multipoles 
/ < 2 are aligned with an axis normal to their alignment plane being oriented along the preferred frame deter- 
mined by the dipole. The preferred direction is displayed by the velocity field of the bulk flow relative to the 
expansion frame of the universe. The parameters are tightly constrained thus the derived modifications lead 
to robust predictions for testing our theory. The 'dark' flow was recently discovered by Kashlinsky et al. to 
be about 700km/ s which seems in good agreement with our predictions for the induced dipole of order 3pK. 
Placed in this context, the discovery of the bulk flow by Kashlinsky et al. becomes even more interesting as it 
may provide a probe of the preinflationary physics and a window onto the landscape multiverse. 

PACS numbers: 98.80.Qc, 11.25.Wx 
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I. INTRODUCTION 

The birth of our universe is one of the most important mys- 
teries left in our understanding of nature. It has also been 
one of the most notoriously difficult problems to address. The 
program we developed over the last few years lU 0> HI H Htl 
advocates a new approach to this central problem. Our investi- 
gation of the selection of the initial conditions of the universe 
differs radically from the previous attempts. This program is 
based on the point of view that: the initial system must be 
studied as a non-equilibrium event; the quantum dynamics of 
the gravitational degrees of freedom needs to be incorporated; 
and, our physical theories have to be extended to a multiverse 
framework. 

Here we show how the observed 'dark' flow is naturally 
predicted in this theory to be a remnant of the birth of the 
universe from the landscape multiverse. Our theory is based 
on the sole assumption that quantum mechanics is valid at all 
energy scales. The calculation is carried out by proposing to 
allow the wavefunction of the universe to propagate on the 
landscape multiverse [1]. Its cornerstone achievement is that it 
provides a sensible answer to the question: 'why did our uni- 
verse start from such an extremely ordered state?' ||3lltl- The 
dynamics of the wavefunction leads to a superselection rule 
that selects 'survivor' universes to be those born at energies 
high enough, (ordered states), that survive the backreaction 
of matter modes, and 'terminal' universes to be those initial 
states that start at sufficiently low energies that can not survive 
the collapse induced by the backreaction of matter modes onto 
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the initial patch. The working model we used 101 for demon- 
strating the application of our proposal, was the landscape of 
string theory since string theory is the leading candidate for 
quantum gravity at present. 

In order to derive testable predictions for this theory, in 
2006 we embarked upon an investigation of remnants in 
present day observables, of the preinflationary dynamics that 
led to the birth of our universe from the landscape multiverse. 
The handle we used to track down the preinflationary rem- 
nants of the matter backreaction and nonlocal entanglement 
with modes beyond the Hubble volume, was based on the 
unitary evolution of the universe. According to unitarity, in- 
formation about the preinflationary remnants of superhorizon 
sized entanglement and non-equilibrium dynamics of this sys- 
tem on the landscape can not be lost. Such investigation IHQl 
led to a series of predictions for traces on CMB and LSS of 
remnants of the preinflationary dynamics of our state on the 
landsacape. Two of these predictions, namely: a giant void at 
z ~ 1 and size about 200 Mpc, as well as a suppressed am- 
plitude of the CMB fluctuation,induced by the superhorizon 
non-inflationary inhomogeneities, ag — 0.8, were succesfully 
tested recently by radio measurements |6] and, WMAP -5- 
year data |7]. A by-product of this analysis were the predic- 
tions for another void near the quadrupole, and of the align- 
ment of the lowest CMB multipoles. These remnants origi- 
nate from the same mechanism, the 'tilting' of the potential 
induced by the superhorizon entanglement |l3l0,|5|l- 

The scope of this Letter is to show that the recent observa- 
tion of the 'dark' flow liSJ tests the third prediction we made in 
||3l namely, the existence of superhorizon sized non-Gaussian 
and noninflationary inhomogeneities that arise from the back- 
reaction of superhorizon matter modes and nonlocal entangle- 
ment of our universe with all other patches and massive modes 
beyond the horizon. Considering that the characteristic length 
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of the noninflationary inhomogeneities Li is 10^ times larger 
than the horizon r^, then their main contribution is to add 
power to the lowest CMB multipoles. For clarity we explic- 
itly calculated the contribution of these superhorizon inhomo- 
geneities to the quadrupole in |3, 4]. We showed they produce 
an induced quadrupole in the CMB, ^, which is added to the 
inflationary one. The noninflationary modification to the grav- 
itational potential by these superhorizon inhomogeneities was 
also derived, (Eq. (3.4) in i4|] ), and as shown there |U,|1,|3] 
this derivation makes no allowance for phenomenology, since 
its parameters are tightly constrained. We derived the quan- 
tum coherence length Ly |3], a length that describes the dis- 
tance at which the quantum interference effects between our 
Hubble volume and modes beyond the horizon become signif- 
icant and of order unity relative to the inflationary density per- 
turbations. For example, for GUT scale inflation, the expres- 
sion for L\ calculated in ^ gives a coherence length of order 
Li ~ 900r// with rn the horizon radius. The cumulative effect 
of the superhorizon noninflationary inhomogeneities thus has 
a characteristic wavenumber and scale k\ ~ \/L\. Besides the 
quadrupole, the nonlocal modification to the gravitational po- 
tential by these modes also induce a dipole in the CMB. The 
dipole can be estimated within our model of the landscape 
multiverse as shown in Sec. 3. The rest frame determined by 
the noninflationary dipole does not coincide with the expan- 
sion frame, rather it induces relative motion with respect to 
the expansion of the universe. The noninflationary dipole has 
the same amplitude as the one derived for the quadrupole in 
j3] but it is a factor of 0(Li/r//) larger than the quadrupole. 
We estimate the bulk flow arising from the induced dipole pre- 
dicted here 11,0] and compare its value to current data. 

Kashlinsky et al fs*] have recently claimed a stunning dis- 
covery: they observe just such coherent cosmological bulk 
flow of matter aligned along the dipole, on scales out to 
at least 300/i^' Mpc, and with amplitudes of order 600 — 
1000 km/sec by using the kinetic Sunay'ev-Zel'dovich effect 
(KSZ). Such a directional flow relative to the expansion of 
the universe frame, derived by comparison with the cosmic 
X— ray background, provides evidence that the induced dipole 
in the CMB is of noninflationary origin. They see no sign of 
convergence of such a bulk flow which indicates that the co- 
herence length is at least horizon size. Thus the name: 'the 
dark flow'. Since the correlation length of such a flow is hori- 
zon size then the noninflationary perturbation that gives rise to 
the dipole has to be induced by superhorizon inhomogeneities. 
But, considering that the flow provides a new preferred frame 
relative to the expansion then these inhomogeneities can not 
be of inflationary origin. 

These results are extremely difficult to accomodate within 
the paradigm of gravitational instability in the concordance 
ACDM cosmology. In this cosmology, the pecuUar bulk ve- 
locity would decrease as Vims o= r^' on a scale r and should 
be extremely small on the scales probed in Ref.|8]. Besides, 
the observed dipole in an inflationary universe is simply kine- 
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matic, i.e. due to our local relative motion with respect to 
the galactic plane. The intrinsic CMB dipole is canceled out 
by the gravitational potential term, therefore unobservable. In 
other words in concordance cosmology, the rest frame defined 
by the expansion of the universe coincides with the CMB rest 
frame. The directional motion of clusters aligned with the 
dipole and with horizon size correlation length, i.e. the 'dark 
flow' observed by |8], provides a second frame relative to the 
expansion of the universe frame and indicates that corrections 
to the gravitational potential for the induced dipole must be 
due to superhorizon sized noninflationary inhomogeneities. 

The modification to the gravitational potential in the uni- 
verse, produced by remnants of the quantum nonlocal entan- 
glement, and the subsequent quadrupole derived in 101, are 
reviewed in Sec. 2. The corresponding noninflationary dipole 
induced by the 'tilted' potential is discussed in Sec. 3. The 
predicted imprints trace out the non-local entanglement ef- 
fects among different horizon regions and modes beyond our 
horizon, that started in different sectors of the landscape and 
interfere with one another to induce decoherence. The deriva- 
tion of the decoherence effects on the wave function of the 
Universe was done via the usual mechanism of tracing out 
superhorizon fluctuations to construct an effective theory that 
describes inflationary physics within the given horizon. 

Our contention, then, is that these observations of bulk flow 
can be construed as evidence for the birth of the universe from 
the landscape multiverse imprinted on the superhorizon sized 
nonlocal quantum entanglement between our horizon patch 
and others that began from the landscape. When we calculate 
the size of the induced dipole in our theory and convert it into 
a bulk velocity dispersion, we will see that for the constrained 
values of our parameters we arrive at a velocity dispersion of 
order ^ 670 km/sec, remarkably close to the observed value 
of -700 km/sec ! 



II. REVIEW: BIRTH OF THE UNIVERSE FROM THE 
LANDSCAPE MULTIVERSE 

We will only sketch out the model we use here so as to set 
the stage for our explanation of the dipole anisotropy and refer 
the reader to Refs.[2, 5, 9] for more details. 

Think of the panoply of all possible ground states for the 
true theory of everything as a lattice of potential wells. We 
proposed in [ 1 ] to let the wavefunction of the universe to prop- 
agate on the landscape of these potential wells. The super- 
space for the wavefunction is determined by the values of the 
various moduli fields needed to specify the vacua, as well as 
on the gravitational variables which is the scale factor a[t] of 
the 3— geometries, LI 01 : these variables form the configura- 
tion space for the wavefunctions. If we speciaUze to the string 
landscape, we can separate the lattice sites into those that de- 
scribe supersymmetric vacua and those in which supersym- 
metry (SUSY) is broken. The wavefunction follows a Bloch 
wave pattern on the SUSY sector thus it is not of interest to 
cosmology since decoherence is hard to achieve. But in the 
broken SUSY sector, the vacuum energies are randomly dis- 
tributed, as shown by Douglas and Denef |9] within a certain 
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range given by the string scale. We treat it as a disordered 
lattice. The strength of the disorder is of order the string scale 
or the Planck scale, since the fundamental parameter of the 
theory determines the range of vacua energy distribution. 

Disordered lattices have the interesting possibility of An- 
derson localization 1 11 ] as proposed in . This allows for the 
localization of the wavefunction around each one of the bro- 
ken SUSY lattice sites, which in turn allows us to think of each 
lattice site as the potential starting point of new "baby uni- 
verse" in the sense that different horizon regions could have 
started from different vacua characterized by the physics of 
the given vacuum site. For this reason, we proposed in Ollsll 
to view the landscape of string theory as the physical phase 
space for the initial states and view the ensemble of all wave- 
function solutions born out of this landscape as the physical 
multiverse predicted by the landscape of string theory. How- 
ever, since the full wave function of the Universe contains all 
that information, the physics in any given site must be entan- 
gled with that of all other lattice sites. Each initial wavepacket 
in the family of solutions for the wavefunction of the universe 
contains matter and vacuum energy. Decoherence among such 
solutions is induced by the backreaction of massive perturba- 
tion around the vacua site where the wavepacket is localized. 
With decoherence taken into account, the quantum dynami- 
cal evolution of matter and gravitational degrees of freedom 
of the wavepacket leads to the following superselection rule: 
the initial states that contain a vacuum energy high enough to 
overcome collapse from the backreaction of matter modes re- 
sult in 'survivor' universes, whereby the quantum to classical 
transition for this universe is succesful. Such initial domains 
can grow to physically large classical universes. But initial 
wavepackets that start at low energies can not overcome the 
collapse induced by the backreaction of matter and thus lead 
to 'terminal' universes since these states do not survive the 
transition to the classical phase and can not grow. The latter 
are wiped out of the phase space of initial conditions, i.e. from 
the ensemble of the initial states which construes the multi- 
verse. One of the 'survivor' universes leads to the birth of our 
universe, The subset of the original ensemble selected as 
'survivor' universes comprises the present multiverse. 

The initial states that led to the birth of physical universes 
follow a unitary evolution. It is known that a mixed state can 
not evolve into a pure state under a unitary evolution |3, 4], a 
fact that ensures that traces of the initial entaglement among 
different patches and of the backreaction of superhorizon mat- 
ter modes onto our patch, do survive in the present observable 
sky. The cumulative effect of the superhorizon nonlocal en- 
tanglement among the various domains and fields has a num- 
ber of physical consequences. First, in order to understand 
the physics inside a given horizon, superhorizon fluctuations 
present in the wave function of the universe must be traced 
out, generating a reduced density matrix for the relevant vari- 
ables. As usual in a reduced density matrix, the off-diagonal 
coherences measure how entangled the system is, and in par- 
ticular, they give a measure of the scale L\ on which these 
cumulative entanglement effects are large. Another way to un- 
derstand this scale, denoted by L\ is to think of it as the mini- 
mum distance we need to go far out of the horizon upon which 



quantum interference effects with other domains and super- 
horizon modes become sufficiently significant for revealing 
the quantum nature of the universe. Equivalently, L\ is the 
scale at which quantum interference corrections becomes of 
the same order as the classical curvature inhomogeneities. 

We showed in llH HI] how the tracing out of superhorizon 
modes gives rise to corrections to the Friedmann equation for 
the scalar field that drives inflation. These corrections are not 
arbitrary, but are in fact related to parameters specifying the 
broken supersymmetry sector of the landscape. We have com- 
puted these corrections in previous work, and list them below. 
The modified Friedmann equation for the inflaton field with a 
potential V{^) on our domain, being embedded in the multi- 
verse, now becomes: 
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Here we have taken %iiGn = M^^, = and b is the 

SUSY breaking scale, which in our picture also corresponds 
to the width of the Anderson localized wavepackets centered 
about a vacua site of the landscape. We can understand the 
nonlocal modification term, F{b,V), to the Friedmann equa- 
tion, by noting that the entanglement and backreaction of 
other patches and modes on our initial state induces a shift 
in the energy of the wavepacket of our universe, correspnd- 
ing to a deviation from its classical path in the phase space. 
This information is contained in the phase of the perturbed 
wavepacket, (see Eq.(3.3-3.8) in |3] for the derivation of this 
modification). 

As mentioned above, the off-diagonal elements of the re- 
duced density matrix give us information about the amount of 
entanglement present between vacua. In particular, the char- 
acteristic length scale Li of the quantum interference, which 
contains the information about entanglement with all the dif- 
ferent domains and superhorison modes, was derived in [3] in 
terms of the parameters of the model: 
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As can be seen by these expressions and their derivation in 
there is no room allowed for phenomenology. 



III. 'DARK' FLOW, SUPERHORIZON ENTANGLEMENT 
AND THE INDUCED DIPOLE 

We reviewed how tracable remnants for testing our theory 
for the birth of the universe, are induced by the superhorizon 
sized inhomogeneities with a characteristis scale given by Li 
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that arise from the nonlocal entanglement of our domain with 
modes and patches beyond the horizon. These noninflation- 
ary superhorizon inhomogeneities are added to the inflation- 
ary CMB, leading to a modification of the Newtonian potential 
<I> Eq l3.1l derived in [3, 4]. The nonlocal modification to the 
gravitational potential of the universe is relevant here because 
the 'tilt' induced by the modification in the potential drives the 
bulk flow as we now demonstrate. The 'tilted' potential with 
the subsequent residual (noninflationary) quadrupole induced 
by it, were calculated in Ref. ISii^J to be: 



Quadrupole Temperature in \iK 



(3.1) 



where <I> is the inflationary, standard gravitational poten- 
tial before modifications 5<I> are included, and the induced 
quadrupole anisotropy resulting from this potential. 
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From the 'tilt' modification term 5<t> in the potential <I>, 
Eqs |3.1| 13.21 it is straighforward to derive the expression for 
the induced dipole. 
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The nonlocal effects only modify the largest scales. We can 
see this from Eq. ( 13.11 ). since the ratio r/Li encoded in 5<I>, 
with L\ ^ r//, can significantly affect power only to the CMB 
multipoles of the largest scales, i.e. for r ~ rn = //^'. In fact 
as can be seen from our dipole expression, even the ampli- 
tude of the induced quadrupole Z = 2 is already suppressed 
by a factor of order (9(r///Li) as compared to the dipole. 
This is equivalent to the fact that since the distance scale for 
our noninflationary perturbations is Li ~ l/ki then k-x I. 
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Thus, since 5p exp(— ;^ ■ x) ~ 5p[l — ik ■ x — ^^'2 4 
multipole decomposition of the density perturbations will be 
dominated by the lowest multipoles. We perceive such a long 
wavelength perturbation as a linear gradient field because the 
induced multipoles are dominated by the first correction term, 
k-x i.e. the induced dipole, with the next correction being the 
quadrupole ^ (k-x)^. Unlike the null primordial CMB dipole 
produced by inflationary perturbations, the lack of cancelation 
between the gradient term (first term) and the 'potential' term 
(second term) in the Sachs Wolfe equation Eq |4.1| results in 
the nonzero dipole 1 12] of Eq l3.3l Thus the predicted dipole 
and its subsequent bulk flow provide a preferred frame. Since 
most of the power goes to the dipole, then naturally the align- 
ment plane of the lowest induced multipoles, < 2, ^ (k ■ x)' , 
is along the gradient of the induced dipole frame. 

The natural alignment of the low multipoles for / < 2, with 
their axes normal to the frame determined by the dominant 
dipole (k-x), proves useful for discriminating our predictions 
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FIG. 1 : The bulk velocities in km/s generated by entanglement ef- 
fects and the temperature of the induced quadrupole and dipole in 

: 3.810"''Mp/ and 



/jK for representative values of the parameters b : 
Vo of order the GUT scale. 



by making use of polarization maps. The predicted dipole in 
Eq l3.3l provides a preferred frame and it gives rise to a bulk 
flow with velocity Eq l3.4l The flow along this dipole pre- 
ferred frame is relative to the rest frame of the expansion of 
the universe. The correlation length of such bulk flow, Eq |3.4| 
induced by the noninflationary dipole is of order as we 
show below. We can now readily estimate the peculiar ve- 
locity field from the induced dipole temperature anisotropics 
through the simple relation 
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where the calibration factor a is related to the averaged 
optical depth of clusters T. We use a calibration factor of 
IQQkm/s for each QA/jK anisotropy in converting the dipole 
temperature of Eq. l3.3l to a peculiar velocity for the bulk flow 
Eq l3.4l This flow has horizon size correlation. Our predicted 
values for the bulk flow (Eq l3.4l l. the dipole (Eq l3.3b and, the 
quadrupole (Eq |3.2| ). are shown in Fig.l. For GUT scale in- 
flation Li is of order 900///. As can be seen from Fig.l, the 
induced dipole temperature predicted here is about 3/j/r, and 
the induced quadrupole is suppressed by the factor H^^ /Li 
thus it is about 900— times less than the dipole. 
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FIG. 2: Likelihood Contours constraining 'b' as a function of V{(S^) 
for values of primordial spectrum /\ at k = 0.002 Mpc are shown. 
The value of X is fixed as A. = 0. 1 in this figure. 



It is worth emphasizing that the derived expressions for the 
induced quadrupole and dipole above can not be assigned dif- 
ferent values by a change of parameters. The parameter 'b' 
corresponding to the initial width of the wavepacket (SUSY 
breaking scale), is tightly constrained by an upper and lower 
bound relative to the inflation energy scale V ((])), derived in 
lilt]. The lower bound is illustrated by the contour plots of 
Fig. 2. These bounds constrain the ratio of the parameters /j/V, 
that enter the expression for the quadrupole and octupole, to 
be a constant, thereby implying that the predictions made by 
this theory are robust and can be used to test it. 

We show the results of our derivation for the bulk flow 
driven by the induced dipole in Fig.l. Amazingly, the bulk 
flow predicted by this theory is of order 650 — 750^m/i, a 
value in agreement with the range observed by |8], (the exact 
value depends on the calibration factor used) and the predicted 
dipole is of order 3/j/r. The calibration used in Fig.l corre- 
sponds to IQQkm/s flow for a dipole temperature OA/jK. In 
our theory the induced dipole is of noninflationary origin thus 
its bulk flow provides a preferred direction aligned with it, that 



moves relative to the expansion frame. This result is in agree- 
ment with the 'dark' flow measurements of [8J, deduced by 
comparing the bulk flow field with an all-sky cosmic X-ray 
map which is at rest with the expansion frame. 

An interesting question related to the robustness of our pre- 
diction is: would our prediction for the dipole change if we 
use a different set of values for {b,V), within their allowed 
range? As can be seen from the likelihood plot of Fig. 2 which 
uses our current limits for the density perturbations to be less 
than 10"^, ( Pk ~ 10"'°, the blue line in Fig.2), the values 
of b are constrained by the inflaton energy V to be such that 
the ratio b/V remains a constant. The 'plateau' in Fig.2 is 
excluded by the upper bounds on 'b' derived from the slow- 
roll conditions in |3]. But since all the nonlocal modifica- 
tions terms involve the ratio of these two scales b/V instead 
of their values, then these correction terms predict the same 
strength independent of the chosen values for b,V. (It is sim- 
ple to understand why all the correction terms induced to: the 
Friedmann equation /^[/JjV]; to the Newtonian gravitational 
potential 5<I> Eq l3.1l to the cumulative coherence length Li in 
Eq l2.2l and for the induced dipole Eq. 13.31 and quadrupole 
Eq l3.2l are always given as function of the ratio of the two 
parameters b and V. These correction terms are nothing more 
than 'greybody' factors on the DeSitter temperature of the in- 
flationary domain given by the inflation energy scale. Since 
in our theory we have two fundamental scale, 'V' and 'b', then 
that will give rise to 'greybody' correction terms to the zero 
order (nearly) DeSitter temperature 'V'. Greybody corrections 
are always a function of the ratio of the two scales since they 
compare the influence of the 'bath' (the landscape multiverse 
for our case) onto the system (our Hubble volume) thus the 
reason why modifications are functions of b/V). We had not 
until now appreciated in [1, 4J the significance of our previ- 
ous constraints for the constant ratio of b/V. To conclude, 
a constant ratio (see Fig.2 and f3\) between these scales in 
Eqs |3. 11331 implies that our predicted values for the induced 
dipole, quadrupole and bulk flow, remain the same indepen- 
dently of the energy scale of inflation V in our domain. This 
conclusion leads to a robust prediction for the 'dark flow' and 
the dipole since their values are centered around those shown 
in Fig.l. 

Moreover in this framework due to the alignment of the 
induced low multipoles £ < I, the octupole-quadrupole axes 
normal to their alignment plane, has to be along the same 
direction as the bulk flow determined by the dominant mul- 
tipole, the dipole. This leads to the very interesting conse- 
quence that the preferred axis addressed in 11 H, (an ob- 
servation discussed by many authors lITsll ). is not actually a 
preferred axial symmetry but instead it is a preferred direc- 
tion defined by the induced dipole. The preferred frame, as 
demonstrated by the relative flow (the 'dark' flow) of all clus- 
ters, including our own local group, with respect to the rest 
frame defined by the expansion of the universe, is determined 
by the simple fact that our Hubble domain is embbedded in 
a much larger environment, the landscape multiverse. The 
existence of a preferred frame can have major consequences 
in our interpretation of current data. One such consequence 
is related to the Doppler shifted photons, due to our 'dark' 
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flow, whereby we would perceive a north-south asymmetry 
with one side appearing 'hot' and the other 'cold'. The char- 
acteristic scale encoding the information about the cumula- 
tive nonlocal entanglement with everything beyond the hori- 
zon, determines the characteristic wavelength of the induced 
superhorizon inhomogeneities Li . The superhorizon inhomo- 
geneities and entanglement encoded on CMB and large scale 
structure (LSS) described here and in |3, 4] can account for 
the list of observed anomalies. Through the predictions made 
in our theory, which include: the lack of power at the largest 
scales in temperature anisotropics; the alignement of the low- 
est multipoles £ < 2 to the dipole gradient; the existence of a 
preferred frame given by the induced dark flow of the dipole 
relative to the universe's expansion rest frame; the north-south 
asymmetry; and, the unusual high power on the correlation 
with polarization Cf^ IHtl, we might be glimpsing at a 
more coherent picture of cosmology in a multiverse frame- 
work. These testable predictions could place the WMAPI?*] 
and NASAliSJ findings into context. 



IV. IS THE DARK FLOW A SIGNAL FROM THE 
LANDSCAPE MULTIVERSE? 

Before we discuss tests for discriminating the predictions 
made by our theory from other mechanisms, let us first briefly 
review whether other models offer plausible explanation for 
the bulk flow observed by [SJ. 

The concordance cosmology could in principle produce 
a bulk flow via unusual superhorizon curvature inhomo- 
geneities, as shown in |12]. Besides, the problem with Vrms 
mentioned in Sec.l, there are a few more issues with this ap- 
proach: first, such superhorizon fluctuations have to be tuned 
to be less than the density perturbations of (9(10^^) but there 
is no physical framework or motivation for producing pertur- 
bations that are specially tuned to have a specific value, but no 
less and no more; second, the more serious problem with such 
a scenario lies in their detection. Even if we accidentally hap- 
pened to have inflationary superhorizon perturbations, of just 
the right amount to produce a IQQkm/s flow, we would still 
not be able to know of its existence because the spatial hy- 
persurface defined by the isotropy of CMB coincides in this 
case with the isotropy of expansion, as shown in lfl3l . The 
latter implies that there is no relative flow of one frame with 
respect to the other. Mathematically the inflationary dipole is 
canceled by the dipole anistropy of the background Newto- 
nian potential as follows: according to the Sachs-Wolfe effect 
lfl2i I13I1 the temperature anisotropics in the direction n by an 
observer today at comoving coordinate r = are: 



Ar(r = 0;n) 



+ 



(4.1) 



Here R, E denote the reception and emission events respec- 
tively. The first term could in principle give rise to a dipole 



anisotropy due to the presence of the gradient operator. How- 
ever, in the standard inflationary picture the Newtonian poten- 
tial <I> is directly related to the induced quadrupole, thus the 
result of Ref.| 12] that the inflationary induced dipole is can- 
celed between the first and second term. Furthermore, a co- 
hesive picture of cosmology that is in agreement with present 
data, requires that besides an accidental 'tilting' of the grav- 
itational potential <I> by just the right ammount for the flow, 
we also have: an 'accidental' alignment of the lowest multi- 
poles along a preferred direction; an anomalously high power 
on the polarization spectra; but an anomalously low power for 
the temperature spectra at the largest scales. So many acci- 
dental alignments make the likelihood of agreement between 
a concordance cosmology that doesn't allow for noninflation- 
ary perturbations, and the best fit to data, extremely improba- 
ble & 12, M, i5J • 

This is the essential difference between the standard sce- 
nario with our approach; the gravitational potential <I> gains 
a contribution from a highly nontrivial second source which 
is independent of the inflationary curvature perturbations and 
is a remnant of the superhorizon size nonlocal entanglement 
from the preinflationary dynamics of our domain in the land- 
scape. For this reason, the induced dipole Eq l3. 31 predicted 
by our theory of the Initial Conditions, has a different ampli- 
tude from the one produced by the second term which ensures 
that the cancelation does not occur. As we showed, the expan- 
sion frame and the CMB frame do not coincide for this case. 
Such a net result for the dipole and the quadrupole produces a 
bulk flow relative to the expansion of the universe frame. The 
induced quadrupole of Eqn l4.1l can be used to discriminate 
among models, as we briefly describe below. 

Another possible explanation for the observed dark flow 
could be provided by scenarios of isocurvature fluctuations 
llT2ll or other models with scalar field perturbations ITtIi . 
Isocurvature fluctuations arise from fluctuations of other 
scalar fields, such as axions or some other PNGB's. Other- 
wise, the growth and the dynamics of subhorizon and super- 
horizon isocurvature perturbations is just like that of any other 
scalar field, including that of inflaton perturbations. Since 
isocurvature fluctuations are independent of curvature fluctu- 
ations, the cancelation between the two dipole contributions 
does not occur in Eq l4.1l for the same reason as the lack of 
cancelation in our theory. But, in order to get the bulk flow 
with horizon size correlation among velocities of clusters, su- 
perhorizon size fluctuations are needed. Isocurvature modes 
are highly constrained by current data 1 7], to be such that they 
do not contribute to the energy perturbations in the early uni- 
verse. Some degree of unnaturalness is required in order to 
minimize or exclude the subhorizon isocurvature modes and 
demand that from the onset of inflation these modes were al- 
ways superhorizon and remain so today. One could imagine 
there may be a scenario where the subhorizon modes are cut- 
toff and the tight data constraints are fulfilled by some special 
isocurvature model. The question thus becomes: can such 
models give rise to the 'tilted' universe? Well, surprisingly, 
due to the expansion of the universe being dominated by dark 
energy the answer to this question seems to be no for the fol- 
lowing reason: since isocurvature fluctuations are constrained 
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not to contribute to density perturbations, then these fluctua- 
tions can grow and excite the background radiation field, only 
during an axion-dominated universe. In complete analogy 
with the growth of inflaton perturbations, isocurvature modes 
are linearly proportional to the scale factor a[t] and can grow 
only during a matter dominated universe,. Therefore although 
their contribution is constrained to be negligible for most of 
the life of the universe, these modes could still have chance 
to grow and dominate at late times if the universe becomes 
axion-matter dominated. If that were the case, isocurvature 
modes would grow and dominate as a result of which a com- 
pensating temperature anisotropy would arise at the dipole 
level, 8T /T\ciip ~ —5^/3. Such a dipole could, in principle, 
produce the observed dark flow [8]. As we know the universe, 
is not axion dominated but rather dark energy dominated. The 
problems is that the growth of perturbations stops during the 
dark energy domination era. Therefore the isocurvature per- 
turbations can not grow and dominate the energy density of 
the universe. The condition for exciting the radiation field in 
order to produce the noninflationary dipole is not fulfilled in a 
dark energy dominated universe. 

An important issue which deserves further investigation, re- 
mains: how can we discriminate our predictions here from 
those of other models? Let us briefly sketch out the follow- 
ing test with polarization maps. We will present this method 
in detail elsewhere |18]. Polarization of CMB is generated 
through rescattering of the temperature quadrupole. Most of 
it, Eprim, is generated at the last scattering surface, with a sec- 
ondary contribution, Eis^, owed to the evolution of the poten- 
tial <I> via the ISW effect. For this reason, polarization maps 
have proven useful for probing dark energy via the ISW ef- 
fect encoded in Eisw [19]. But, when a bulk motion with 
peculiar velocity v relative to the rest frame of expansion is 
present, then there is also a kinematic component of polar- 
ization, Ek, since a quadrupole, proportional to v^, is being 
generated. The polarization distribution traces out the under- 
lying quadrupole distribution. Thus, for the components Eis^ 
and Eprim which trace out the primordial quadrupole produced 
by inflaton or another scalar, we should expect a random dis- 
tribution in the polarization skymap. But in our theory, the 



induced quadrupole proportional to is oriented with re- 
spect to a preferred frame determined by the dipole and the 
transverse velocity field of clusters v. Since the correlation in 
the bulk flow is horizon size then the kinematic component 
of the polarization map Ek that traces out this orientation of 
the quadrupole will not be random. If we can extrapolate Ek 
from the skymaps we are able to test this theory. This can be 
done by substracting Eisw and Eprim from the total measured 
polarization, which removes their random distribution and al- 
lows us to deduce whether Ek is random or oriented. Should 
this mechanism be the right one for addressing the 'dark' flow 
then the kinetic component of the polarization skymap will 
not have a random distribution 11 8.1 . Related to this test is the 
fact that the spectral intensity of the kinetic polarization map 
will be frequency dependent thus provide further information 
on the predicted alignment of the quadrupole and dipole gradi- 
ent. Once we have better probes for the Ek components, these 
tests will also be relevant for providing more accurate probes 
of dark energy via the ISW effect imprinted on the polariza- 
tion Eisw Upcoming data from Planck mission will increase 
the precision of polarization maps and perhaps allow the ex- 
trapolation of the kinetic component. 

What is surprising is that with these parameter choices, 
we arrive at bulk flow velocities consistent with recent 
observations! 8] and coherent with all the observed anomalies. 
Should this theory ultimately be tested to put the findings of [? 
] and |7] in context, then such discovery is providing us with a 
sneak preview of the vast landscape beyond our horizon and, 
a Copernican Principle extended to the whole universe 12011 . 
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